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Abstract 

It  is  shown  by  surface  lattice  dynamics  that  a  new  type  of  surface  vibrational 
resonance  arises  in  those  frequency  regions  where  there  is  a  strong  depletion 
in  the  bulk  phonon  density  of  states.  The  presence  of  these  pseudoband  gaps 
is  due  to  the  higher  Fourier  components  in  the  phonon  dispersion  relations 
introduced  by  the  particular  coordination  of  atoms  in  layers  parallel  to  the 
surface.  A  quantitative  analysis  based  on  surface  lattice  dynamics  of  the 
recorded  electron  energy  loss  spectra  of  Cu  and  Ni  suggests  that  the  outer¬ 
most  surface  interlayer  force  constant  is  rather  close  to  the  bulk  value.  This 
resonance  is  found  to  exist  throughout  the  T.f  direction  and  makes  an 
avoided  crossing  with  a  resonance  derived  from  a  band  gap  at  the  i’-point. 
An  explanation  is  given  for  the  disappearance  of  the  divergent  van  Hove 
singularities  in  the  projected  bulk  density  of  states  upon  projection  on  a 
surface  layer. 


1.  Introduction 

The  vibrational  properties  of  clean  metal  surfaces  have 
recently  attracted  a  lot  of  attention  both  from  an  experimen¬ 
tal  and  theoretical  point  ot  view.  These  studies  have  been 
made  possible  by  new  surface  sensitive  vibrational  spec¬ 
troscopies  such  as  inelastic  He  scattering  and  electron  energy 
loss  spectroscopy  (EELS).  These  techniques  probe  vibrations 
with  such  high  frequencies  that  the  surface  lattice  dynamics 
cannot  be  described  fully  within  the  continuum  elasticity 
theory.  In  this  frequency  regime  the  atoms  vibrate  with  large 
relative  displacements  such  that  the  surface  geometric  struc¬ 
ture  and  the  surface  interatomic  forces  are  expected  to  play 
a  decisive  role.  Recent  measurements  on  Ni(100)  by  EELS 
[1,2]  and  Ag(l  1 1)  by  inelastic  He  scattering  [3,  4]  demon¬ 
strated  that  such  information  can  be  extracted  from  surface 
vibrational  spectra. 

The  possibility  to  observe  dipole  active  surface  vibrational 
modes  on  clean  metal  surfaces  was  recently  demonstrated  for 
the  (100)  surfaces  of  Cu  and  Ni  by  EELS  [5].  In  a  recent  letter 
we  reported  the  observation  of  dipole  active  resonance  modes 
on  the  (110)  surfaces  of  Cu  and  Ni  [6].  These  modes  were 
shown  from  surface  lattice  dynamics  calculations  to  be  a  new 
kind  of  resonance  arising  from  a  pseudoband  gap  in  the 
density  of  states  for  longitudinal  phonons  propagating 
normal  to  the  surface.  This  gap  defines  a  region  where  the 
bulk  phonon  density  of  states  is  significantly  depleted  and  has 
a  simple  structural  explanation  in  terms  of  the  coordination 
of  the  atoms  in  the  layers  normal  to  the  (1 10)  direction  of  fee 
crystals.  A  recent  analysis  of  the  bcc(l  1 1)  surface  has  shown 
that  this  type  of  surface  resonance  is  a  general  effect  which 
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results  from  higher  harmonic  terms  in  the  bulk  phonon  dis¬ 
persion  relations  introduced  by  the  particular  coordination  of 
atoms  in  the  layers  parallel  to  the  surface  [7], 

In  this  paper  we  show  in  detail  how  a  surface  vibrational 
resonance  arises  in  situations  with  pseudoband  gaps  in  the 
bulk  phonon  density  of  states.  The  surface  is  observed  to  split 
off  a  mode  from  the  region  of  high  density  of  bulk  phonon 
states  into  the  pseudoband  gap  region  where  the  strong 
depletion  of  the  density  of  states  causes  the  mode  to  become 
a  resonance.  An  argument  is  also  given  for  explaining  why  the 
divergent  van  Hove  singularities  which  are  present  in  the 
projected  bulk  density  of  states  disappear  for  a  projection  on 
the  surface  layers.  These  divergences  disappear  because  an 
incident  phonon  at  those  frequencies  interfere  destructively 
with  the  phonons  scattered  from  the  surface.  On  the  [110} 
surfaces  of  fee  Ni  and  Cu  crystals  the  pseudoband  gap  is 
shown  to  exist  along  the  redirection  in  the  surface  Brillouin 
zone  (SBZ).  The  corresponding  surface  vibrational  resonance 
makes  an  avoided  crossing  with  a  resonance  derived  from  a 
gap  mode  at  the  e-point  in  the  SBZ. 

A  detailed  quantitative  comparison  of  the  measured 
spectra  with  the  calculated  EEL  spectra  show  that  the  loss 
peak  derived  from  the  resonance  is  reproduced  with  a  value  for 
the  outermost  surface  interlayer  force  constant  within  ±  1 5% 
from  its  bulk  value.  In  contrast,  large  oscillatory  relaxations 
of  the  surface  layers  were  observed  by  low  energy  electron 
diffraction  (LEED)  for  these  surfaces  [8],  The  EEL  spectra 
also  give  information  about  the  dipole  activity  of  the  surface 
layers.  The  shape  of  the  loss  spectrum  is  well  reproduced  by 
only  the  two  outermost  surface  layers  giving  the  dominant 
contribution  to  the  dipole  activity.  The  dipole  strength  is 
found  to  be  of  the  same  order  of  magnitude  as  measured  for 
the  (100)  surfaces  [5], 

The  vibrational  structures  of  the  surface  appear  in  the 
measured  spectra  through  a  specific  projection  of  the  surface 
vibrational  density  of  states.  While  the  relevant  projection 
probed  in  inelastic  He  scattering  is  essentially  given  by  the 
displacements  of  the  outermost  surface  atoms  normal  to  the 
surface  [9],  the  relative  rigid  displacements  of  the  surface 
layers  is  the  relevant  projection  in  inelastic  dipole  scattering 
[5].  Such  vibrational  density  of  states  have  been  evaluated 
here  from  surface  lattice  dynamics  for  a  semi-infinite  sub¬ 
strate  using  simple  force  constant  models. 

The  force  constant  models  have  been  extracted  from  Born 
von  Karman  analysis  of  inelastic  neutron  scattering  data  with 
particular  emphasis  on  the  high  frequency  vibrations.  The 
bulk  phonon  dispersions  of  Cu  and  Ni  are  well  described  by 
a  single  nearest  neighbouring  force  constant  [10.  IIJ.  The 
surface  force  constants  have  been  chosen  to  have  the  same 
values  as  the  bulk  layers.  This  choioce  should  be  viewed  as 


yy.yv; 


Geometric  Structure,  Pseudoband  Gaps,  and  Surface  Vibrational  Resonances  on  Metal  Surfaces 


a  reference  model  with  no  more  justification  than  that  it 
describes  the  large  effect  of  the  loss  of  coordination  of  atoms 
at  the  surface. 

The  surface  vibrational  density  of  states  has  been  calcu¬ 
lated  using  a  Green’s  function  technique  proposed  by  Lee 
and  Joannopoulos  [12].  This  technique  is  based  on  the 
transfer  matrix  and  its  application  to  surface  lattice  dynamics 
is  described  in  the  Appendix. 

2.  Experimental  procedures 

The  experiments  were  performed  in  a  multitechnique  ultra- 
high  vacuum  system  which  is  evacuated  by  turbomolecular. 
ion.  and  titanium  sublimation  pumps  with  a  base  pressure  of 
4  x  10  "  Torr  [13],  The  electron  spectrometer  is  based  on  a 
double  pass  127°  cylindrical  electrostatic  deflector  for  the 
monochromator  and  analyzer.  The  scattering  geometry  is 
fixed  with  a  total  scattering  angle  of  120°.  The  angular  accept¬ 
ance  of  the  analyzer  is  1.8°  at  full-width-half-maximum 
(FWHM).  The  scattering  plane  containing  the  incident  and 
scattered  electrons  is  defined  by  the  surface  normal  and  the 
[lTO]  crystallographic  direction  for  the  Cu(l  10)  and  Ni(l  10) 
surfaces.  The  spectra  were  recorded  in  the  specular  direction 
at  a  temperature  of  300  K.  Impact  energies  of  3.2  and  4.3  eV 
were  used  for  Cu  and  Ni  respectively. 

The  samples,  which  were  approximately  1  cm  diam  and 
I  mm  thick  disks,  were  cleaned  by  neon  ion  bombardment 
(500eV)  and  annealing  to  750  and  1050  K  for  Cu  and  Ni, 
respectively.  The  samples  were  spotwelded  to  a  manipulator 
with  a  pair  of  0.5  mm  Ta  wires  for  the  Ni  samples  and  with 
0.5  mm  Pt  wires  for  the  Cu  sample.  The  clean  surfaces  dis¬ 
played  sharp  1  x  l  LEED  patterns  without  any  sign  of 
typical  impurity  vibrations,  such  as  O,  C,  or  S,  in  the  electron 
energy  loss  spectra. 

The  vibrational  spectra  of  the  clean  Cu  and  Ni(l  10)  sur¬ 
faces  are  shown  in  Fig.  1 .  A  single  sharp  vibrational  loss  peak 
is  observed  at  20  and  24meV  on  the  Cu(110)  and  Ni(l  10) 
surfaces,  respectively.  Off-specular  measurements  show  these 
losses  to  be  excited  by  the  dipole  scattering  mechanism  [14], 
Energy  gain  peaks  are  also  observed  with  an  intensity  ratio  to 
the  energy  loss  peaks  determined  by  the  Bose -Einstein  distri¬ 
bution  factor  at  300  K.  An  important  feature  of  these  losses 
is  that  the  ratio  of  their  energies  scales  as  the  ratio  of  the 
maximum  bulk  phonon  frequency  which  is  29.7  and  36.7  meV 
for  Cu  [10]  and  Ni  [11],  respectively.  This  fact  suggests  that 
these  losses  are  derived  from  longitudinal  phonons  propagat¬ 
ing  normal  to  the  surface. 

3.  Pseudoband  gaps  and  surface  resonances 

The  notion  of  a  pseudoband  gap  is  illustrated  by  the  disper¬ 
sion  of  longitudinal  bulk  phonons  in  the  [110]  direction. 
A  detailed  analysis  of  the  displacement  fields  for  phonons 
scattered  from  the  surface  shows  how  a  surface  vibrational 
resonance  can  develop  in  this  situation  and  why  the  divergent 
van  Hove  singularities  present  in  the  projected  bulk  phonon 
density  of  states  disappear  on  the  surface  projection.  The 
development  of  a  resonance  on  the  (1 10)  surface  is  contrasted 
with  the  (100)  and  (III)  surfaces  where  no  such  resonances 
appear  within  the  nearest  neighbor  central  force  constant 
model. 


Energy  Loss  (meV) 

Fig.  I.  Electron  energy  loss  spectra  of  the  clean  ( 1 10)  surfaces  of  Cu  and  Ni, 
The  specra  were  recorded  in  the  specular  direction  at  300  K  The  sharp  peak) 
observed  at  20  and  24  meV.  for  Cu  and  Ni  respectively,  correspond  to  surface 


vibrational  resonances. 


the  layers  of  atoms  normal  to  the  surface  can  be  dipole  active. 
Such  a  motion  of  the  bulk  layers  corresponds  to  longitudinal 
bulk  phonons  propagating  in  the  [1 10]  direction.  It  is  known 
since  the  early  studies  by  inelastic  neutron  scattering  that  the 
full  bulk  phonon  dispersions  for  Cu  and  Ni  can  be  well 
described  by  a  Born-von  Karman  model  of  lattice  dynamics 
based  on  central  nearest  neighboring  force  constants  [10,  1 1], 
In  this  model  the  eigenvalue  problem  for  the  longitudinal 
bulk  phonons  propagating  in  the  [110]  direction  is  given  by 


m2wL  =  —  (4wt  -  -  W[  _ i  -  h',  ,2  -  w,  ,). 

and  is  in  the  [100]  and  [111]  directions  given  by 


<o:wL  =  ^(2w, ,  -  w,*,  -  wL  ,), 


where  wL  is  the  displacement  of  an  atom  in  the  Lth  layer  in 
a  direction  normal  to  the  layer,  to  the  frequency,  x  the  central 
nearest  neighboring  force  constant,  and  M  the  mass  of  an 
atom. 

From  the  translational  symmetry  of  the  bulk  layers  the 
solutions  to  eqs.  ( I )  and  (2)  are  simple  plane  waves  u  ,  =  e'''L 
with  energies  uKC)  satisfying  the  disperston  relation  lor  the 
[110]  direction. 


>‘(C)  =  y  j^sin2  j  +  sin:(*0  . 


For  the  (1 10)  surface  of  fee  metals  only  rigid  motions  of  and  for  the  [  1 00]  and  [III]  directions. 
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where  t%  =  8r  .V/.  The  resulting  phonon  dispersion  from 
this  simple  force  constant  model  for  the  lattice  dynamics 
agrees  fairly  well  with  data  from  inelastic  neutron  scattering 
for  all  the  major  crystallographic  directions  as  shown  in 
Fig.  2  when  «„  is  adjusted  so  that  uAi  =  I)  is  equal  to 
the  measured  value  in  the  [110]  direction,  i.e.,  a»(y  =  I)  = 
v  2.  ui0  =  29.7  and  36.7  meV  for  Cu  and  Ni,  respectively. 
The  tit  is  particularly  good  in  the  high  energy  region  and  the 
largest  discrepancies  are  typically  found  in,  for  this  analysis, 
the  less  interesting  low  energy  region. 

A  characteristic  feature  of  the  dispersion  in  the  [1 10]  direc¬ 
tion  is  its  non-monotonic  behavior  with  a  maximum  at 
=  ,'m.  well  within  the  Brillouin  zone.  This  behavior  can  be 
understood  simply  from  the  coordination  of  atoms  in  this 
direction  shown  in  Fig.  3.  In  this  case  an  atom  has  nearest 
neighbors  not  only  in  the  nearest  layer  but  also  in  the  next 
nearest  layer.  The  next  nearest  neighboring  atom  lies  in  the 
[110]  direction  and  causes  the  next  nearest  interlayer  force 
constant  to  be  as  strong  as  the  nearest  interlayer  force  con- 
sum  in  eq.  ( I ).  This  causes  for  instance  the  restoring  force  for 
the  displacement  fields  at  =  0  5  to  be  stronger  than  at 
=  I  0.  In  the  other  two  crystallographic  directions  an  atom 
has  only  nearest  neighboring  atoms  in  the  nearest  neighbor¬ 
ing  layer  and  this  causes  the  dispersion  to  be  monotonic  with 
>  as  evidenced  by  eqs.  (2)  and  (4). 

In  surface  vibrational  spectroscopy  one  probes  the  vibra- 


Ft*  .1  Slruclureof  fee  crystals  in  the  [I  I0|  and  (100)  directions  (a)  Top  view 
i>l  atoms  of  the  (110)  surface  together  with  the  crystallographic  directions 
The  1 1 10)  surface  Brillouin  zone  is  depicted  in  (b)  The  coordination  of  atoms 
in  the  bulk  layers  normal  to  the  surface  are  shown  for  (c)  the  (1 10)  and  Id) 
the  (I00|  directions. 


tionai  density  of  slates  in  the  surface  region.  Before  present¬ 
ing  results  for  the  surface  density  of  states  we  will  consider 
first  the  density  of  states  for  a  bulk  layer.  The  phonon  density 
of  state  g(ui)  for  a  bulk  layer  is  simply  determined  by 

!  dr.o  I  1 

g<W)  =  —  •  <5) 

I  I ; . 

where  the  sum  is  over  all  positive  C,  satisfying  oj  =  t.o(^). 
From  the  phonon  dispersion  relations  given  by  eqs.  (3)  and 
(4)  tf(ro)  can  be  evaluated  straightforwardly  and  is  presented 
for  Ni  in  Fig.  4  (upper  panel).  In  the  low  energy  limit  ui  = 
(,C-t  d and  g(<>)  goes  to  a  constant.  g(a>)  ^  dic^n,  where  c ,  is 
the  longitudinal  sound  velocity  and  d  the  interlayer  spacing. 
In  one-dimensional  problems,  as  in  the  case  considered  here, 
the  stationary  points  in  the  dispersion  relation  defined  by 
d w  dC  =  0  give  rise  to  divergent  van  Hove  singularities  in 
.?(";)[  1 5].  The  divergences  are  in  most  cases  power  singularities 
with  an  exponent  -  I  2.  This  kind  of  singular  behavior  is 
readily  shown  from  the  fact  that  it  is  possible  to  make  a 
Taylor  expansion  o( O  ~  <■>  +  y(;  -  2  around  an 

isolated  stationary  point  C  *  C.  and  eq.  (5)  gives  directly  that 

1  ,,  . 

if(«u)  i - a)  —  to  I  ",  s,  ■  (6) 

s  -I7i 

However,  in  some  exceptional  cases,  which  are  not  encoun¬ 
tered  here.  /  =  0  and  the  next  leading  term  in  the  Taylor 
expansion  gives  rise  to  another  value  for  the  exponent.  For 
metals  ii  is  also  possible  to  have  non-analytical  behavior, 
Kohn  anomalies,  from  the  long-range  interactions  intro¬ 
duced  by  I  he  sharp  Fermi  surface  Tlu-sc  singularities  are  not 
discernible  for  Cu  and  Ni.  The  most  important  point  to  be 
made  here  about  g(<o)  is  the  fact  that  the  [110]  dispersion 
relation  has  a  relatively  large  density  of  states  in  a  rather 
narrow  region  in  to.  22  <  h  $  33  meV.  compared  to  the  low 
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Fig  4  The  density  of  states  *(<o)  for  longitudinal  phonons  projected  on  bulk 
and  surface  layers  for  Ni.  The  results  for  g(<ul  in  the  [1 10],  (111)  and  (100) 
crystal  directions  have  been  calculated  using  the  same  force  constant  model 
as  in  Fig  l 


energy  region  0  <  h  <  22meV.  This  latter  region  will  for 
that  reason  be  called  a  psuedoband  gap.  No  such  region  can 
be  defined  for  the  other  two  directions  [110]  and  [111]. 

In  order  to  evaluate  g(«)  for  a  surface  layer  one  needs  a 
model  for  the  effects  of  the  surface  on  the  force  constants.  An 
obvious  effect  of  forming  a  surface  is  the  corresponding  loss 
of  coordination  of  the  surface  atoms.  Here  we  will  only 
account  for  that  effect  on  the  interlayer  force  constant.  The 
resulting  equations  describing  the  surface  layers  in  the  [110] 
direction  are  given  by. 

,  X 

"r'tt'i  =  —  ( 2h-|  -  w,  -  wy) 

M 

(?) 

r  , 

<•>'  u\  =  —  (3w,  -  w,  -  tv,  -  tv4), 

V/ 

and  tor  the  outermost  surface  layers  in  the  [100],  [III] 
direction 


Note  that  within  this  formalism  the  relaxation  of  the  surface 
layers  is  taken  into  effect  only  through  modifications  of  the 
surface  interlayer  force  constants. 

The  presence  of  the  surface  breaks  the  translational  sym¬ 
metry  and  the  solutions  to  eqs.  (I)  and  (7)  can  no  longer  be 
written  as  a  single  plane  wave.  Rather  a  wave  e  '"'L  incident 
on  the  surface  will  be  reflected  and  can  couple  to  another 
wave  e*'**1  with  the  same  energy  oi(y)  =  This  kind  of 
argument  suggests  the  following  ansatz  for  the  scattered  wave 
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in  the  [1 10]  direction, 

wL  =  e  ,<L  +  R(0  c‘<L  +  R(Q  e’"!t.  (9) 

This  form  for  the  ansatz  is  justified  in  the  Appendix.  Due 
to  the  non-monotonic  behavior  of  the  dispersion  some 
care  is  needed  to  get  the  right  boundary  conditions.  The 
reduced  wavevector  (  has  to  be  chosen  from  the  ranges 
-1  $  y  <  -£mand0  ^  C  <  where  the  maximum  in  the 
phonon  dispersion  occurs  at  ±Cm  In  these  ranges  the  group 
velocity  is  positive,  dui/d(  >  0,  so  that  e  ~"'L  and  e";i  are 
incoming  and  outgoing  waves,  respectively.  The  other  reflected 
wave  e+*,K’L  js  propagating  for  oi/ai0  ^  1/^/2  and  similarly 
has  to  be  chosen  from  the  ranges  where  the  group  velocity  is 
positive,  dtu/d£  >  0.  For  smaller  energies  <  1  / v  2,  ^  is 
complex  and  e'“'L  is  an  evanscent  wave  and  the  sign  of  Im  ’ 
has  to  be  chosen  so  that  it  is  a  decaying  wave. 

The  two  reflection  coefficients  R(0  and  R(y)  are  now 
determined  from  the  two  equations  of  motion  for  the  surface 
layers,  eq.  (7).  The  two  equations  for  R( 0  and  R(£)  are  given 
by. 

d,(0R(0  +  </,(’)£(;)  =  1-1.2  do) 

where 

4(0  =  (v:(C)  -  1.4)  +  ‘(e";  +  e‘:,;) 
and 

4(0  =  (v’(C)  -  3/8)  +  >(e  +  e‘";  +  e,:";). 

Here  v  is  the  reduced  energy  v  =  <oiid„.  The  surface  vibra¬ 
tional  density  of  states  defined  as. 

g(a>)  =  |o'  dC  -  a»(0]  (II) 

can  now  be  evaluated  from  eq.  (9).  It  turns  out,  however,  that 
it  is  more  elegant  and  practical  to  evaluate  g(w)  by  the 
transfer  matrix  method  described  in  detail  in  the  Appendix. 
The  resulting  g(a>)  calculated  by  the  transfer  matrix  method 
for  the  [1 10]  direction  is  depicted  in  Fig.  4  (lower  panel). 

A  noteworthy  feature  of  g( to)  is  that  the  divergent  van 
Hove  singularities  have  disappeared  in  the  projection  on  the 
outermost  surface  layer.  This  can  be  shown  rather  easily  from 
eqs.  (9)  and  ( 10)  to  be  due  to  the  fact  that  one  gets  destructive 
interference  between  incident  and  reflected  waves  resulting  in 
tv,  s  Oat  the  stationary  points.  For  instance  at  s'  =  the 
ansatz  degenerates  to 

UY  =  [I  +  «(Cm)|e  +  «(Cm)e‘";™'.  (12) 

Instead  of  having  two  inhomogeneous  equations  for  K(2)  and 
from  eq.  ( 10)  we  have  now  two  homogeneous  equations 
fori  +  R(sm)  and  /?(£).  These  two  equations  will  in  general 
have  a  trivial  solution  except  in  those  accidental  cases  where 
the  determinant  is  identically  zero.  Thus  the  contribution  of 
g(w)from£  =  Cmis£(tu)  =  |iv,  (£)|J  ,jdw/d£|.  From  eq.  ( 10)  it 
is  evident  that  H’,(0  is  analytic  around  C  =  C™.und  »,(,')  i 
/.(2  -  Cm)  since  iv,(Cm)  =  0.  Similarly  for  ru(C).  w<0  -  <■■>«,  4- 

VtC  -  s'm):/2,  and 


,?M  - 


-  w)11. 


Thus  the  divergent  van  Hove  singularity  ~(<o  -  <■>„,)  '  at 

a  bulk  layer  for  the  one-dimensional  model  turns  into  a 
bounded  van  Hove  singularity  ~(w  --  <um)  '  on  a  surface 


layer.  This  argument  indicates  also  that  the  divergent  van 
Hose  singularity  should  not  exist  in  g(w)  for  a  projection  on 
any  layer  for  the  semi-infinite  substrate.  A  closer  analysis 
reveals,  however,  that  the  bulk  density  of  states  are  recovered 
in  the  limit  L  -*  ao.  For  instance,  for  a  layer  far  inside, 
Z.  S>  1,  |wj2  =S  4  sin2[rr(£  -  £J£]  when  £  -*  £m.  Thusg(to) 
will  rise  rapidly  when  going  away  from  tum  and  have  a  large 
maximum  8L/|y|  at  £  -  £ra  =  1/2 L,  arbitrarily  close  to  cum. 
This  argument  for  the  disappearance  of  the  divergent  van 
Hove  singularities  in  a  projection  of  g(o>)  on  a  surface  layer 
can  be  shown  to  apply  to  more  general  situations.  For 
instance,  it  is  not  necessary  that  surface  force  constants  are 
the  same  as  in  the  bulk  region. 

Most  importantly,  g(a>)  shows  a  sharp  narrow  feature 
around  23  mV  as  seen  in  Fig.  4  (lower  panel).  This  feature 
is  now  shown  to  be  a  surface  vibrational  resonance.  In  a 
situation  where  there  is  an  absolute  band  gap  it  is  well  known 
that  the  surface  can  introduce  a  localized  state  split  off  from 
the  band.  In  the  present  case  there  is  no  absolute  band  gap 
rather  a  pseudoband  gap.  The  surface  can  possibly  split  off  a 
state  from  the  band  which  turns  into  a  resonance  by  over¬ 
lapping  with  the  low  density  of  bulk  states  in  the  pseudoband 
gap.  This  expectation  is  confirmed  from  an  analysis  of  the 
reflection  coefficient  for  cj :oj0  <  1  v2.  *  =  1  +  i*  where 
cosh  7t k  =  j[l  +  v  25  -  32(ai/(u0):]  and  the  complex  part 
gives  rise  to  an  evanescent  wave  ( —  1  )L  e  'L .  The  correspond¬ 
ing  reflection  coefficient  £(£)  is  found  to  have  a  simple  pole 
for  complex  v  at  =  0.663  +  i0.047  (for  Ni,  = 

(24.3  •+-  il.73)meV).  The  existence  of  such  a  pole  with  an 
imaginary  part  cu,  relatively  close  to  the  real  axis  justifies 
calling  this  rather  sharp  peak  a  surface  vibrational  resonance. 
Note  that  the  peak  is  quite  asymmetrical  due  to  interference 
with  bulk  states  in  the  depicted  projection,  a  feature  which  is 
typical  for  Fano-resonances  [16], 

For  the  [100]  and  [111]  directions  the  ansatz  for  the 
solution  to  eqs.  (2)  and  (8)  has  a  more  simple  form 

wL  =  &’:L  +  fi(Oe";t  (14) 

for  0  ^  ^  I.  This  ansatz  inserted  into  eq.  (8)  for  the 

surface  layer  gives  a  simple  form  for  the  reflection  coefficient 
£(£)  =  e",; .  As  a  function  of  v  this  reflection  coefficient, 
£(£)  =  1  —  2 v2  -t-  2iv/v:  -  v4  (when  0  <  £  <  I).  has  no 
poles  associated  with  any  resonances.  The  phonon  density  of 
states  g(cu)  projected  on  a  surface  layer  can  now  be  evaluated 
directly  from  eq.  (II)  and  is  given  by 


g(w)  =  -  V71  -  (W/l  >o)2- 


(15) 


This  density  of  states  shows  accordingly  no  surface  vibrational 
resonances  as  depicted  in  the  lower  panel  of  Fig.  4. 


4.  Comparison  with  experiment 

An  attractive  feature  of  EELS  is  the  possibility  to  analyze 
quantitatively  the  measured  dipole  active  losses  [17],  The 
dipole  loss  function  for  longitudinal  bulk  phonons  is  calcu¬ 
lated  for  Ni  and  compared  with  the  measured  spectrum.  The 
sensitivity  of  the  calculated  spectra  to  changes  in  the  surface 
force  constant  and  the  distribution  of  the  dipole  activity 
among  the  surface  layers  are  also  investigated. 

In  a  recent  letter  it  was  shown  both  experimentally  and 
theoretically  that  the  displacements  of  the  outer  layers  of 


metal  atoms  can  give  rise  to  a  long  range  dipole  field  due  to 
incomplete  screening  by  the  conduction  electrons  of  the  elec¬ 
tric  field  from  the  displaced  ion  cores  [5].  The  strength  of  the 
dipole  field  is  described  by  effective  charges  e*  which  relate 
the  normal  component  of  the  dynamic  dipole  moment  pi.  to 
the  rigid  displacements  wL  of  layer  L  normal  to  the  surface 
through, 

M.-  =  T.e>L-  (16) 

L 


Here  we  use  the  same  model  for  ef  as  in  Ref.  [5], 

c*  =  — cf  =  c*  and  c*  =  0,  L  >  2.  (17) 


Note  that  a  rigid  displacement  of  the  metal  atoms  normal  to 
the  surface  cannot  give  rise  to  a  dipole  moment,  i.e., 
c?  =  0.  The  projection  of  the  phonon  density  of  states 
relevant  for  dipole  losses  is  accordingly  given  by 


IX  «**7.(£) 

I  L 


P[tu  -  aj(0J 


(18) 


where  nf  =  with  e*„  =  (Et  ef2)1  \  is  the  normalized 

field  of  effective  charges.  The  spectral  function  5(a>)  for  the 
dipole -dipole  correlation  function  appearing  in  the  energy 
loss  function  is  related  to  g( cu)  through  [18], 


S(w)  =  [I  +  /»(«)]  e*J 


h 

2Mio 


SM. 


(19) 


where  the  mass  M  of  a  metal  atom  appears  in  the  root-mean- 
square  amplitude  h/2Ma>  for  phonons  with  energy  twu  and 
n(jj)  is  the  Bose-Einstein  distribution  factor. 

From  the  inelastic  dipole  scattering  theory  the  inelastic 
current  /,  (to)  of  electrons  collected  in  the  detector  around  the 
specular  direction  after  experiencing  an  energy  loss  ftw  is 
given  with  sufficient  accuracy  by  [17] 

/,(<»)//,„,  =  ,  -”?e - /(£„,  w,  x)S(w).  (20) 

n~AE„  cos  % 

where  /,ol  is  the  total  integrated  intensity  of  the  elastic  peak  in 
the  energy  loss  spectrum,  m  the  electron  mass.  A  the  area  of 
the  surface  primitive  cell,  and  £„  the  kinetic  energy  of  the 
electron  metdent  with  an  angle  i  from  the  surface  normal. 
The  function  /(£n,  a>,  a)  is  given  by  [17] 

/(£„,  a»,  a)  =  (sin:a  -  2  cos2 a )Y 

A-  (sin2a  +  2  cos2a)  In  X.  (21) 

where  Y  =  0];(0]  +  02),  X  =  (0J  +  0]),0l  0„  =  fuo  2£n 
gives  the  angular  extension  of  the  dipole  lobe,  and  0,  the 
half-angle  of  the  detector  aperture.  The  loss  function  depicted 
in  Fig.  5  is  now  obtained  from  eq.  (20)  by  calculating  the 
projected  phonon  density  of  states  defined  in  eq.  (18)  by  the 
transfer  matrix  method  for  the  distribution  of  effective 
charges  given  in  eq.  (17).  The  parameters  a,  0,.  and  £„  are 
determined  from  the  experimental  conditions  described  in 
Section  2.  and  the  experimental  resolution  was  introduced  by 
a  4meV  Gaussian  broadening  of  g(io).  The  total  effective 
charge  e*,  had  to  be  chosen  to  be  0.034e  and  0.039c  for  Cu 
and  Ni,  respectively,  in  order  to  reproduce  the  measured  loss 
intensities  at  300  K..  These  values  are  of  the  same  order  of 
magnitude  as  for  the  value  determined  previously  for  the 
Cu(l00)  surface  [5].  Because  the  resonance  gives  rise  to  a 
rather  sharp  loss  peak  there  has  been  no  particular  need 
to  have  a  detailed  analysis  of  the  contribution  from  the 
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Fix -5.  Caiculated  electron  energy  loss  function  /,(cu)//0  for  Ni(IIO)  at 
300  K  /0  Is  the  maximum  intensity  of  the  elastic  peak  In  the  energy  loss 
spectrum.  Only  the  two  outermost  surface  layers  are  assumed  to  be  dipole 
active  and  the  total  effective  charge  ?*,  has  been  adjusted  to  0.039c  (e  is  the 
free  electron  charge)  in  order  to  reproduce  the  measured  loss  in  Fig.  I.  The 
instrumental  resolution  has  been  introduced  by  a  Gaussian  broadening  of 
4  meV 

electron-hole  pair  excitations  to  subtract  the  background 

[191. 

The  calculated  position  of  24.5  meV  and  the  peak  width  of 
6meV  for  the  surface  vibrational  resonance  are  in  good 
agreement  with  the  measured  values  for  Ni.  Note  that  the 
value  for  this  peak  position  is  about  1  meV  higher  than  for  the 
peak  position  deduced  from  g(a>)  for  the  projection  on  the 
outermost  layer.  This  difference  is  due  to  the  fact  that  the  low 
energy  bulk  phonons  contribute  much  less  to  this  dipole 
active  projection,  which  suppresses  the  asymmetry  of  the 
peak.  The  peak  position  is  thus  closer  to  the  value  for  the  real 
part  of  the  pole  in  the  complex  co-plane  as  given  in  the 
previous  section.  For  Cu,  g(co)  is  obtained  in  the  advocated 
force  constant  model  simply  by  scaling  the  phonon  energies 
with  a)0(Cu)/w0(Ni)  %  0.81.  This  gives  an  energy  of 
19.8  meV,  in  good  agreement  with  the  measured  value  of 
20meV  observed  in  Fig.  I. 

There  are  no  reasons  to  expect  that  the  only  effect  of  the 
surface  is  the  loss  of  coordination  of  atoms  in  the  surface 
region  as  described  by  eq.  (17).  For  instance,  both  model 
calculations  for  the  total  energy  [20]  and  LEED  measure¬ 
ments  [8]  have  shown  that  the  atoms  relax  oscillatorily  in  the 
surface  region  for  many  metals.  In  these  new  equilibrium 
positions  for  the  atoms  the  force  constants  can  be  different 
from  the  bulk  values.  Off-specular  measurements  of  the 
Rayleigh  surface  phonon  dispersion  on  Ni(l00)  by  EELS 
have  suggested  that  the  interlayer  force  constant  between  the 
first  and  second  layer  is  about  20%  larger  than  the  bulk  value 
[I.  2).  On  Ag(l  1 1)  the  observation  by  inelastic  He  scattering 
of  a  surface  vibrational  resonance  away  from  the  T-point  in 
the  SBZ  could  be  accounted  for  by  a  reduction  of  about  50% 


of  the  radial  part  of  the  force  constant  between  atoms  in  the 
surface  region  [3,  4]. 

The  position  of  the  surface  vibrational  resonance  in  the 
observed  energy  loss  spectrum  should  also  contain  infor¬ 
mation  about  the  surface  force  constants.  The  sensitivity  of 
the  position  of  the  resonance  to  changes  in  the  surface  inter¬ 
layer  force  constants  has  been  investigated  by  calculating  a 
dipole  active  projection  of  the  phonon  density  of  states  g(w) 
for  different  values  of  the  surface  interlayer  force  constants 
for  Ni.  The  resonance  is  predominantly  localized  in  the  two 
outermost  surface  layers  and  should  accordingly  be  most 
sensitive  to  changes  of  the  interlayer  force  constants  within 
these  layers.  The  corresponding  modification  of  the  interlayer 
force  constants  is  described  by  the  following  equation  of 
motion  for  the  surface  layers, 

.V/ai2  w,  =  a,(w,  -  w,)  +  x2(w,  -  w3) 

Ma>1w2  =  3|(w,  —  w,)  +  a(2w,  -  vv,  -  w4)  (22) 

A/arw,  =  a2(wj  —  iv, )  +  a(3wj  —  w,  —  —  w,) 

The  resulting  position  and  width  of  the  resonance  are  tabu¬ 
lated  in  Table  I  for  a  considerable  range  of  values  for  R ,  = 

a,  /a  and  R2  =  a. ,/«.  No  attempt  has  been  made  to  relate 

these  parameters  to  the  observed  oscillatory  relaxations  of 
this  surface  since  that  would  require  a  detailed  microscopic 
model  for  the  interaction  potential.  It  would  be  possible  to 
proceed  with  the  simple  pair  potential,  but  the  results 
obtained  would  be  highly  doubtful  for  a  transition  metal  due 
to  the  different  characters  of  the  contributions  to  the  interac¬ 
tion  energy  from  the  free  electron  like  s-electrons  and  the 
tightly  bound  ^-electrons.  That  the  position  of  the  resonance 
is  much  more  dependent  on  /?,  than  R2  is  no  surprise  since  the 
displacement  field  of  the  resonance  has  its  dominant  weight 
on  the  relative  displacement  of  the  two  outermost  surface 

Table  I.  The  influence  of  the  surface  interlayer  force  constants 
on  the  resonance  position  and  width.  The  surface  interlayer 
force  constants  /?,  =  a, /a  and  R2  =  a2/a  are  normalized  to 
the  hulk  interlayer  force  constant  i.  The  resonance  position  w„., 
and  the  width  Tr„  are  normalized  to  the  position 
(  =  24.5  meV  for  Ni)  and  width  I"?„  (  =  3.1  me  V  for  Ni.  esti¬ 
mated  full  width  at  half  maximum)  corresponding  to  the  situa¬ 
tion  a,  =  a2  =  a.  The  area  enclosed  by  the  solid  line  is  the 
range  of  values  in  accordance  with  the  observed  resonance 
position  when  the  error  bars  are  taken  into  account.  The  dashed 
entries  for  the  resonance  widths  indicate  that  the  resonance  is 
no  longer  well  defined 
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Users.  The  experimental  resolution  is  such  that  it  can  deter¬ 
mine  the  peak  position  within  I  meV.  The  corresponding 
range  tor  the  position  of  the  observed  resonance  relative  to 
the  calculated  position  is  then  0.94  to  1 .02  and  is  enclosed  by 
a  solid  line  in  Table  I.  A  large  range  of  values  for  R:  is 
acceptable  but  not  for  /?,.  The  Table  suggests  that  /?,  lies 
between  0.8  and  1.1.  There  are  also  notable  changes  in  the 
w  idth  of  the  resonance  when  Rt  and  R2  are  changed.  For  large 
values  of  R,  the  resonance  emerges  into  the  region  of  large 
bulk  phonon  density  of  states  and  the  width  is  no  longer  well 
defined.  With  the  present  resolution  of  electron  energy  loss 
spectrometers  it  is  not  meaningful  to  extract  any  information 
about  these  parameters  from  the  observed  resonance  width. 

Another  consideration  to  be  taken  into  account  is  how  the 
effective  charges  are  distributed  among  the  surface  layers.  In 
the  case  of  the  Cu(IOO)  surface  the  results  from  a  jellium 
model  calculation  for  Cu  suggested  the  distribution  defined 
in  eq.  (17)  [5].  The  application  of  the  same  model  for  the 
Cull  10)  surface  gives,  however,  that  even  the  third  and 
fourth  lasers  have  an  appreciable  effective  charge  ef  mainly 
due  to  a  smaller  interlayer  distance  in  this  direction.  In  Fig.  6 
(lower  panel)  we  present  results  for  g(oi)  calculated  for  Ni 
with  two  different  distributions  for  =  ef  e*,  extending  to 
the  third  and  fourth  layers  and  compare  with  the  result  from 
the  distribution  defined  in  eq.  (17).  For  the  other  two  distri¬ 
butions  the  resonance  peak  is  still  prominent  and  does  not 
change  its  position,  but  the  strength  of  the  states  in  the 
upper  bulk  band  region  has  been  appreciably  enhanced.  The 
measured  loss  spectra  for  Cu  and  Ni  do  not  indicate  such  a 
strong  contribution  from  the  bulk  states. 

Thus  our  analysis  of  the  loss  spectra  suggests  that  there  are 
no  dramatic  changes  in  the  surface  interlayer  force  constants 
from  the  bulk  values. 

5.  Dispersion  of  the  resonance  along  the  T-.V  direction 

The  dispersion  of  surface  vibrational  modes  along  different 
directions  in  the  SBZ  has  been  shown  to  be  feasible  to 
measure  for  a  few  metal  surfaces  by  inelastic  He  scattering  [3! 
and  off-specular  EELS  [I],  Therefore  it  is  of  interest  to  know 
how  the  resonance  disperses  away  from  the  T-point.  It  is 
found  that  the  resonance  exists  and  is  derived  from  a  pseudo- 
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band  gap  even  out  to  the  .V-point.  The  resonance  makes  an 
avoided  crossing  with  another  resonance  derived  from  a 
surface  phonon  in  a  bulk  band  gap.  Close  to  the  T-pomt  the 
resonance  leaves  the  bulk  subbands  and  appears  as  a  surface 
phonon. 

Along  the  r~jf  direction  the  displacements  of  the  atoms 
partition  into  two  classes  due  to  the  reflection  plane  sym¬ 
metry.  The  odd  modes  are  polarized  in  the  v-direction  and  are 
symmetry  forbidden  to  couple  with  displacements  of  atoms 
polarized  in  the  x-z  plane  which  form  the  even  class.  In  the 
nearest  neighboring  central  force  constant  model  the  motion 
of  atoms  in  the  v-direction  gives  rise  to  a  monotonic  dis¬ 
persion  of  the  corresponding  phonons  with  no  pseudoband 
gaps.  Henceforth  the  .v-motion  will  not  be  considered  further 
The  equations  of  motion  for  displacements  of  atoms  in  the 
v-direction  and  in  the  e-direction  are  coupled  along  the  T-.V 
direction  and  are  for  the  bulk  layers  given  by 


uru,  =  -[(4 


2  cos  ( nc,))uL  -  cos  (it;  2)(u,  ,  +  u, . , ) 


+  i  sin  (^/2)(wL.,  -  wL  ,)], 

,  1 

w*»’t  =  77  t4lv'.  ~  cos(Kf/2)<HY*,  +  wL-\) 


(23) 


-  (wt,;  +  wL  J  +  i  sin  (nc  2)(ul.i  -  u,  ,)]. 

where  u,  e'4'*’  and  wL  e'4'*’  are  displacements  in  the  v- 
and  e-directions,  respectively,  of  an  atom  at  position  R  in 
Layer  L.  and  ;  =  k^ai^ln  is  the  reduced  wavevector  along 
the  direction.  At  the  point  r-point(;  =  0)  the  equations 
of  motion  for  uL  and  \vL  are  decoupled  and  eq.  ( I )  is  recovered 
for  tv,. 

The  solutions  to  eq.  (23)  are  plane  waves  u,  =  u[z.  7)  e  ;/ 
and  tv,  =  w(c,  7)  c'<L  which  would  result  in  two  branches  of 
the  dispersion  relation,  one  lower  <nt  =  (:>,(;.  7)  and  one 
upper  w„  =  <;>„(<;,  7).  The  behavior  of  these  two  branches  at 
7  =  0.6  is  illustrated  in  Fig.  7.  If  one  artificially  removes  the 
coupling  between  u,  and  tv,  then  the  dispersion  for  phonons 
polarized  in  the  v-direction  is  monotonic  and  crosses  twice 
the  dispersion  for  phonons  polarized  in  the  r-direction.  The 
latter  dispersion  is  non-monotonic  due  to  the  strong  coupling 
to  the  second  nearest  neighboring  layer.  The  coupling  present 
in  eq.  (23)  between  u,  and  w,  causes  these  two  branches  to 
make  two  avoided  crossings  with  corresponding  interchange 
of  character  and  makes  them  both  non-monotonic  w  ith  7 

The  influence  of  the  surface  on  the  force  constants  is 
modeled  in  the  same  way  as  in  Section  3  by  taking  into 
account  only  the  loss  of  coordination  of  atoms  in  the  surface 
region.  In  this  complex  case  we  will  not  attempt  to  write 
out  the  form  of  the  scattered  wave  for  an  incident  wave  It 
is  much  more  tractable  to  generate  results  for  the  surface 
vibrational  density  of  states  by  using  the  transfer  matrix 
method.  This  method  cannot  be  applied  directly  to  this  system, 
however,  due  to  the  fact  that  the  dynamical  submatrix  D 
between  the  principal  layers  is  singular  (see  Appendix)  This 
matrix  D,„  can  be  regularized,  however,  by  introducing  a 
small  second  layer  coupling  (x  \t)rtu. .  .  +  m,  •)  into  the 
equations  of  motion  for  «,  in  eq.  (2^1  Ihevaluesolr  =  uni 
is  found  to  be  sufficiently  small  for  an  accurate  calculation  ol 
the  phonon  density  of  states.  This  value  lor  r  is  much  smaller 
than  the  errors  in  the  nearest  neighboring  force  constant 
model  used  to  describe  Cu  and  Ni 
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F'e  '  dispersion  of  bulk  phonons  in  the  |llll|  direction  of  Ni  and  the 
corresponding  surface  phonon  density  ol  states  eli  <|  at  ;  =  k,a  =  0  6 

The  left  panel  shows  the  two  branches,  an  upper  "U"  and  a  lower  "L" 
branch,  ot  the  dispersion  in  the  (l  10|  direction  tsolid  lines)  arising  from  an 
avoided  crossing  between  phonons  polarized  in  the  x ■  and  .--directions, 
respectively.  The  dashed  lines  show  the  dispersion  when  the  interaction 
between  these  two  polarizations  is  turned  off  The  right  panel  shows  the 
phonon  density  ol  states  gout  protected  on  the  .--motion  (solid  line)  and  the 
i -motion  (dashed  line)  of  the  outermost  surface  layer 

The  results  in  Fig.  7  for  the  phonon  density  of  states 
g((.d.  ;)at;  =  0.6  projected  on  the  x-motion  and  the  r-motion 
ot  the  outermost  layer  show  several  prominent  features. 
There  is  a  localized  state  in  the  gap  at  1 1.1  meV.  the  surface 
phonon  S,  in  the  notation  of  Ret.  [2 1  ],  being  split  off  from  the 
hulk  subbands  by  the  reduction  ot  the  restoring  forces  in  the 
surface  region.  The  v-projeetion  of  ?(w.  2)  shows  a  narrow 
peak  at  24  meV  just  below  the  minimum  energy  of  the  upper 
branch  which  can  he  interpreted  as  a  state  being  split  off  from 
the  upper  branch  and  turning  into  a  resonance  due  to  overlap 
with  states  in  the  lower  branch.  Thus  the  origin  of  this 
resonance  is  the  same  as  for  the  resonance  discussed  in  the 
work  on  Ag(  III)  where  an  "anomalous"  peak  was  observed 
in  inelastic  He  scattering  [3],  There  is,  however,  another 
narrow  peak  in  the  r-projection  of  g(w,  2)  around  I9meV. 
The  non-monotonic  behavior  of  the  lower  branch  suEgcsts 
the  interpretation  that  this  peak  is  a  resonance  derived  from 
the  corresponding  pseudoband  gap  of  the  lower  branch 
below  wt  (2,  C  =  I).  The  upper  branch  shows  similar  non¬ 
monotonic  behavior  with  a  pseudoband  gap  in  between  <»„  mi„ 
and  t'zu(2,  2=1)  which  results  in  a  resonance  at  31.1  meV. 
very  close  to  w„(2.  2=1)-  However,  its  dominant  ampli¬ 
tudes  are  on  layers  further  inside  the  surface. 

By  calculating  the  x-  and  r-projections  of  g(</z.  2)  on  the 
outermost  layer  for  several  values  of  2  between  0  and  I  the 
behavior  of  the  surface  vibrational  modes  can  be  followed 
along  the  T  Y direction  as  shown  in  Fig.  X  At  the  ,f-point  we 
have  three  surface  phonons  for  displacements  polarized  in  the 
v  r  plane  (i)  S,  the  Rayleigh  surface  phonon  (u)  S„  which 
exists  only  close  to  .f  und  (in)  S-  a  gap  mode  The  labelling  of 
the  modes  are  taken  from  Ref.  [21]  except  for  S„  which  was  not 
identified  in  their  slab  calculations.  The  mode  S. ,  is  localized  on 
the  second  and  third  layer  and  is  predominantly  polarized  in 
the  .--direction.  This  mode  turns  into  a  resonance  VIS,,  inside 
t  e  lower  bulk  subband  and  lies  in  the  pscudoband  gap  just 
ow  nZ|  (t;,  s  =  I )  At  around  2  =  0.5  0  6  the  resonance 
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h?  X  Dispersion  of  the  resonance  along  the  T.f -direction.  The  dispersion 
ol  the  resonance  arising  from  the  pseudoband  gap  (squares)  makes  an 
avoided  crossing  with  the  dispersion  of  the  resonance  (circles)  derived  from 
the  S-I.V  )  surlace  phonon  The  solid  lines  give  the  maximum  hw,  hm 
and  minimum  hoi,  mm,hto „  energies  of  the  lower  and  the  upper  boundaries 

of  bulk  subbands,  respectively.  The  dashed  lines  show  hoi,  (;,  ;  =  I)  and 
2  =  1 1  for  phonons  with  a  dashed  wavevector  ;  =  k.a  s  2x  where  a 
is  the  lattice  constant.  The  notation  ol  the  surface  phonons  S,  and  S.  are 
taken  from  Ref.  |2I). 

interacts  with  the  MS-  resonance  and  makes  an  avoided 
crossing  with  a  corresponding  interchange  of  character.  The 
resonance  VIS,  is  a  continuation  of  the  gap  mode  S-  into 
the  bulk  subbands  and  becomes  mainly  polarized  in  the  x- 
direction  tor  0.6  <  2  $  1.0.  When  2  approaches  the  T-point 
<2  =  Ol  MS-  goes  over  into  the  resonance  discussed  in  the 
previous  sections  and  is  mainly  polarized  in  the  c-direction. 
from  the  T-poinl  to  the  avoided  crossing  the  width  of  MS- 
remains  roughly  the  same  (about  3.5  meV)  and  after  it  inter¬ 
changes  character  it  sharpens  appreciable  to  a  width  less  than 
0.5  meV .  V1Sn  broadens  and  gets  more  localized  on  the  outer¬ 
most  layer  away  from  the  T-point  and  just  at  the  crossing  the 
width  is  about  2  meV.  After  the  crossing  the  width  remains 
about  the  same  and  sharpens  up  only  just  before  leaving 
the  bulk  subband.  Thus  at  the  crossing  the  widths  of  the 
resonances  overlap,  which  makes  the  avoided  crossing  less 
well  detined 

6.  Summary 

\  new  kind  of  surface  vibrational  resonance  is  shown  from 
surface  lattice  dy  namics  to  exist  on  surfaces  hav me  a  pseud o- 
band  gap  in  the  bulk  phonon  density  ol  states  The  surface 
vplits  oil  a  mode  from  a  region  of  high  densitv  of  states  into 
a  pseudohand  gap  region  where  the  density  of  states  is  lareclv 
depleted  This  behavior  is  illustrated  lor  phonons  having  a 
surlace  component  ol  the  wavevector  along  the  T.V  direction 
in  the  SUZ  of  the  (110)  surfaces  of  Cu  and  Ni  Recent 
measurements  and  analysis  of  the  Fc  hcc  (III)  surface  have 
shown  this  type  of  resonance  to  be  a  general  effect  [7]  In  these 


canes,  the  pseudoband  gap  is  a  geometric  structure  etFect 
caused  by  the  particular  coordination  cf  the  atoms,  which 
leads  to  higher  Fourier  components  in  the  non-monotomc 
bulk  phonon  dispersion  relations. 

At  the  T-point  the  resonance  is  dipole  active  and  has  been 
observed  by  EELS  on  the  ( 1 10)  surface  of  Cu  and  Ni  and  the 
(111)  surface  of  Fe  [7],  From  these  observations  it  has  been 
possible  to  obtain  some  information  on  the  surface  interlayer 
force  constants.  In  particular,  the  positions  of  the  loss  peak 
can  only  be  reproduced  for  Ni  when  the  outermost  surface 
interlayer  force  constant  lies  within  —20%  to  10%  of  the 
bulk  value.  Along  the  f\f  direction  in  the  SBZ  of  Cu  and  Ni 
the  resonance  makes  an  avoided  crossing  with  a  resonance 
derived  from  the  SA-?)  surface  phonon.  This  novel  behavior 
should  be  possible  to  observe  by  inelastic  electron  or  He 
scattering  at  large  parallel  wavevector  transfers. 

Finally,  this  analysis  suggests  in  general  that  this  type  of 
surface  vibrational  resonance  should  be  observable  not  only 
by  inelastic  electron  dipole  scattering  but  by  other  surface 
spectroscopies,  such  as  inelastic  He  scattering,  on  a  variety  of 
surfaces  at  points  in  the  SBZ  where  a  bulk  phonon  dispersion 
is  non-monotonic  and  consequently  has  a  pseudoband  gap. 
The  origin  of  these  effects  is  directly  related  to  the  geometric 
structure  of  the  surface. 
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Appendix 

In  this  Appendix  it  is  shown  how  the  transfer  matrix  method 
proposed  by  Lee  and  Joannopoulos  [12]  can  1  :  applied  to  the 
calculation  of  surface  vibrational  density  of  states.  The 
method  is  illustrated  for  the  surface  lattice  dynamics  problem 
of  longitudinal  phonons  propagating  normal  to  a  fee  (110) 
surface.  Furthermore,  this  method  justifies  the  choice  of  the 
ansatz  for  the  scattered  waves  in  eqs.  (9)  and  (14). 

The  first  step  in  this  method  is  to  form  principal  layers, 

here  labelled  by  an  integer  n.  n  =  1.2 . from  the  layers 

of  atoms  parallel  to  the  surface  such  that  the  dynamical 
matrix  only  introduces  interactions  between  displacement 
fields  in  nearest  neighboring  principal  layers.  In  the  present 
case  two  layers  form  a  principal  layer.  The  column  vector  W„ 
denotes  displacement  fields  in  the  prmcipal  layer  n. 

W,U)  =  :.  i  =  1,2.  (All 

In  terms  of  these  column  vectors  lVn  the  eigenvalue  problem 
for  the  bulk  layers  can  be  written  as, 

.  --  a;,  ik„  ,  =  o,  rr  =  i,2 . 

( A  2 ) 

and  the  corresponding  equation  for  the  surface  layers  is  given 
by. 

<:  -  A,„)fK,  -  A„  IK  =  0.  (A3) 

Here  =  nr  and  A*,,  A,,,  and  A.,,,  are  (2x2)  dynami¬ 
cal  submatrices  formed  from  the  full  dynamical  matrix 
D(L.  1.  )  which  can  be  obtained  directly  from  eqs.  (I)  and 


(7).  Aio(i.  j )  =  A(2L  +  /  -  2.  2L  +  j  -  2),  A, „</.;)  - 
D(2L  +  i  —  2,  2L  +  j).  L  denotes  a  bulk  layer,  and 
A.,,, ( i.  j  )  =  DU.j).  For  instance,  Dm  is  given  by, 


A,  | 


i 

| 


( A4) 


Equation  (A2)  shows  explicitly  that  there  are  only  interac¬ 
tions  between  displacement  fields  in  nearest  neighboring  prin¬ 
cipal  layers.  Since  Am'  exists  W„.{  can  be  directly  expressed 
in  terms  of  the  two  preceding  column  vectors  fK„  and  Wn  ,  by 
a  simple  rearrangement  of  eq.  (A2)  as. 

»K-i  =  A0 , 1  ( ~  -  Dw)Wn  -  An,'  AT,  IK  n  =  2.  3 . 

(A5) 


This  equation  shows  that  it  is  possible  to  construct  a  matrix 
Hr)  which  relates  the  displacement  fields  in  two  principal 
layers  n  +  2andn  ■+■  I  to  the  corresponding  fields  in  the  two 
preceding  principal  layers  n  and  n  —  1. 


IK.; 


( A6) 


The  matrix  T(:)  is  the  transfer  matrix  and  is  given  by  the 
product  of  the  following  two  matrices. 


77--)  = 


A,)i  (c  Ann)  ~  Dn,  Af,  i 


An  I  (c  —  Drn )  —  Adi'A.ii  ' 

I  0 


( A7) 


By  iterating  eq.  (A6).  a  displacement  field  in  any  principal 
layer  can  be  determined  from  their  values  on  the  surface 
layers  as. 


1  \ 

j  IK,  \ 

=  r  <:) 

1 

Equation  (A3)  for  IK,  and  IK,  gives  only  2  equations  for  4 
displacement  fields  and  are  not  sufficient  to  determine  IK,  and 
IV, .  Further  restrictions  are  found  by  introducing  the  appro¬ 
priate  boundary  conditions.  That  can  be  done  by  analyzing 
the  eigenvalues  and  eigenvectors  of  the  dynamical  matrix. 

For  the  bulk  layers  the  solution  of  eq.  (A2)  is  given  by 
translational  symmetry  as  plane  waves. 


( A9) 


where  the  reduced  wavevector  2  satisfies  the  bulk  dispersion 
relation  w  =  <.%[sin‘(rrC  2)  +  sin’(rr2)]  -  as  given  by  eq.  (3) 
in  Section  3.  In  terms  of  the  variable  /.  =  el,;  this  dispersion 
relation  is  equivalent  to  a  polynomial  of  degree  4  in  /.  and  has 
accordingly  4  roots  %.  k  =  1.2,  3.  4.  The  eigenvectors  I[(c) 
of  T(z)  can  now  be  directly  formed  from  this  plane  wave 
solutions. 


KC) 


(A  10) 


and  the  associated  eigenvalue  for  K,(;)  is  given  by  and 
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T(z)Vk  —  Vk .  The  eigenvalues  are  distinct  away  from  the 
critical  points,  dcu/d £  =  0,  and  the  corresponding  eigen¬ 
vectors  span  the  4-dimensional  space  of  displacement  fields  of 
two  adjacent  principal  layers.  Thus  W2  and  Wx  can  then  be 
simultaneously  expanded  in  terms  of  Vk, 
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ckVk. 


(All) 


This  equation  and  eq.  (A8)  give  directly  that  the  displacement 
field  for  any  principal  layer  can  be  expressed  as, 

=  I  cki?Vk,  n  =  0,1,2 -  (A  12) 

k- I 

or  in  terms  of  the  displacement  field  wL  for  a  layer  L, 

=  I  c*  e"'c*i,  (A  13) 

k  »  I 

where  ck  =  This  form  of  solution  in  eq.  (A13) 

justifies  the  ansatz  made  in  eqs.  (9)  and  (14)  in  Section  3.  The 
solution  corresponding  to  scattered  wave  can  be  found  by 
imposing  the  outgoing  boundary  conditions  as  discussed  in 
Section  3.  This  restricts  the  solutions  to  depend  on  two 
parameters.  These  two  parameters  can  then  be  determined 
from  the  two  equations  for  the  surface  layers. 

A  more  convenient  way  to  evaluate  the  vibrational  density 
of  states  g(«,  { nL })  than  using  the  scattered  wave  solutions 
appearing  in  eq.  (A  13)  is  to  determine  first  the  resolvent 
matrix  (a  Green’s  function)  U(L.  L'\  z)  associated  with  the 
dynamical  matrix  D(L,  L').  This  resolvent  is  defined  by, 

£  [:5(L  L")  -  D(L,  L”)]U(L\  L'\  z)  =  S (L,  L'),  (A14) 


and  the  vibrational  density  of  states  is  given  by. 


g(oj.  \nL})  =  --ImY  ntU(L,  L‘\  (to  +  iO  *)2)«t . 

^  L.L 

(A  15) 


The  transfer  matrix  approach  can  now  be  applied  by  con¬ 
sidering  the  resolvent  (2  x  2)  submatrices  (/„„  (r)  with 
respect  to  the  principal  layers  and  they  are  defined  as, 

U„„  ( i,j ■  z)  =  (J(2n  —  L  4-  i,  In'  +  L  -  j\  r). 


(/,  |(c)  simultaneously  in  terms  of  those  eigenvectors  with 
|/t|  <  I.  This  point  and  the  fact  that  for  complex  ;  the 
eigenvectors  are  divided  evenly  into  two  classes  |A*  |  <  I  andi 
UJ  >  1.  respectively,  were  shown  in  detail  for  the  general 
case  in  the  original  work  by  Lee  and  Joannopoulos  [12],  Let 
k  =  I,  2  label  the  two  eigenvectors  with  U»|  <  I  and  intro¬ 
duce  the  two  associated  (2  x  2)  matrices, 


=  Vt(i) 

KUJ)  =  K<‘  +  2). 


(A19) 


The  expansion  of  the  two  submatrices  of  the  resolvent  in 
these  two  eigenvectors  now  becomes, 


=  (WVA\ 

W, .,(=>/  V^L  a) 


(A20) 


where  the  coefficients  in  the  expansion  forms  a  (2  x  2) 
matrix  A.  These  two  resolvents  are  now  specified  by  4 
parameters.  The  4  surface  layer  equations  in  eq.  (A17)  for 
f/n(j)  and  CA,(r)  will  now  completely  determine  these 
parameters.  This  can  be  done  by  first  eliminating  the  matrix 
A  from  eq.  (A20). 


t/,,(r)  =  W^WC'UU(Z).  (A21) 

Furthermore,  by  inserting  this  expression  for  U2i(z)  into 
eq.  (A  17)  a  simple  linear  matrix  equation  is  obtained  for 
(/i  ,(r)  which  can  be  solved  by  a  matrix  inversion, 

Uu(z)  =  (r  -  A*,  -  DnWuWc'y'.  (A22) 

Thus  for  every  frequency  w  the  vibrational  density  of  states 
can  be  evaluated  from  eqs.  (A  15)  and  (A22)  by  diagonaliz- 
ation  of  a  (4  x  4)  complex  matrix  and  by  inversion  of  two 
(2  x  2)  matrices.  £/u(-)  W«H  have  simple  poles  at  those 
frequencies  corresponding  to  localized  vibrational  modes  at 
the  surface.  Similarly,  the  resonances  appear  as  poles  in  the 
complex  frequency  plane  but  not  on  the  physical  Riemann 
sheet  of  C/U(z).  However,  the  other  Riemann  sheets  of  C/, ,  (r) 
should  be  possible  to  construct  from  other  choices  for  the 
eigenvectors  in  eq.  (A  19). 
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i,j  -  L  2.  (A  16) 

To  obtain  the  vibrational  density  of  states  for  the  surface 
layers  it  is  sufficient  to  evaluate  (7U  (z).  The  resolvent  matrix 
element  t/„  ,(r)  satisfies  the  same  equations  as  W„,  eq.  (A2), 
except  at  the  surface  layers  where  the  equations  have  an 
inhomogeneous  term, 

|:-OJI/l.lW-H,t/u(:)  =  1  (A17) 

Similar  to  the  construction  of  W„,  U„  ,(z)  can  be  constructed 
from  L\  ,(c)  and  U1A(z)  by  iterating  the  transfer  matrix. 


U'»,,<--) 


=  T"(z) 


Uu{z) 

U,,(z) 


(A  18) 


Some  care  is  needed  to  get  the  correct  physical  Riemann  sheet 
of  the  resolvent  as  a  function  of  r.  On  this  sheet  U„A(z)  has  to 
be  decaying  with  n  when  imparting  a  small  positive  imaginary 
part  it  to  o»,  z  =  (m  +  ie)2.  Such  a  decay  is  evidently 
achieved  by  expanding  the  two  column  vectors  of  V2A(z)  and 
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